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Abstract: 
We find a new set of Semi-Empirical Mass-Formula coefficients which is a, „ag ,A.,A,, for nuclei with odd mass 
number using least squares fitting method, and we have also added another term for Semi-Empirical Mass 
Formula related to isospin A, using single-particle model. 


When calculating value of standard deviation between values of binding energies of our calculated nuclei and 
experimental values before and after adding isospin term, we found that adding isospin term to Semi-Empirical 
Mass Formula lowers the value of standard deviation. 

And when comparing values of binding energies of our calculated nuclei with values of previous studies, we found 
that our calculated values are better than values calculated in these studies. 
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1. INTRODUCTION 


The liquid drop model is considered one of the first and most important nuclear models, it assumes that the 
nucleus has properties similar to properties of a liquid drop, and by means of this model, the binding energies 
are calculated according to mass number and atomic number by means of the Semi-Empirical Mass Formula 
(SEMF) equation, which is known as the equation of Bethe-Weizsacke[1] . 

According to this equation, nuclear binding energy is a set of terms, each of which expresses a distinct 
characteristic of the nucleus. 

This equation consists of five basic terms: volume energy term, surface energy term, Coulomb energy term, 
asymmetric energy term, and pairing energy term. 

Usually, values of coefficients in SEMF are calculated by proportioning to experimental data using Least Square 
Method, which is considered one of the effective methods for calculations of this type of calculation[2]. 


The aim of this work is to use the method of least square fitting to find a new set of coefficients of the SEMF for 
100 odd-A stable nuclei, that is, they have an odd mass number (odd-even nuclei and even-odd nuclei) using a 
code designed for this purpose by MATLAB, and we adding isospin term using single-particle model. Then we 
compare standard deviation between values of binding energies of our calculated nuclei and experimental 
values before and after adding isospin term. 


2. NUCLEAR BINDING ENERGY 


The sum of the number of protons (which is equal to the atomic number of elements in periodic table) and the 
number of neutrons within the nucleus is called the mass number A= Z +N. 

Protons and neutrons are bound to each other inside a nucleus by very large and short-term force called 
nuclear force, and this force differ from other of known forces in nature such as electrical force and 
gravitational force. As a result of the existence of these large and short-range force, we consider that the 
nucleus is a bound system, so we need an energy of several million electron volts (order of Megaelectron volts 
(MeV)) to separate components of the nucleus compared to several electron volts to separate electrons from 
their atomic orbits (we need an energy of 13.6 eV to separate the hydrogen atom's electron from its orbital) [3]. 


The least amount of energy needed to separate a components of the nucleus of a proton and a neutron is called 
binding energy, and this energy is equal to the energy required to form a nucleus consisting of a proton and a 
neutron separated from each other[3]. Thus, mass of the nucleus is less than sum of masses of its components 


of protons and neutrons, and therefore we can write mass of the nucleus as follows[3] 


AM =ZM , +NM, —M(A,Z) (1) 
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Where My mass of the hydrogen atom, which is equal to 1.007940u, and M, mass of the neutron, which is 


equal to 1.008665u, M is mass of the atom, and therefore binding energy of the nucleus is given as 


E, =|ZM,+NM,-M(A,Z)|c’ (2) 


where c° =931.494061 MeV /u the square of speed of light in a vacuum. 


3. LIQUID DROP MODEL 


The liquid drop model is considered one of the most important nuclear models. This model was proposed by 
Weizacker in 1935 and Bohr developed the basic hypotheses of this model in 1937[1]. 


In this model, the nucleus was likened to a liquid drop, and this analogy was chosen according to the following 
three consideration 


- Nuclear matter is incompressible, so[4] 


V aR? 
1 Vaa (3) 
Rad} 


Where R is a radius of nucleus and A is a number of nucleons, so the volume of nucleus prportion whith 
number of nucleons, as in liquide drop. 
-The approximate stability of binding energy per nucleon for medium and heavy nuclei, and this indicates that 


each nucleon interact with a limit number of nucleons, meaning that the nuclear forces are saturated, as the 
forces acting between the molecules of a liquid. 


-The nuclear force are same for nucleons, which means they are not related to the nucleon being a proton or a 
neutron as in the case of molecules in a liquid drop. 


4. SEMI-EMPIRICAL MASS FORMULA (SEMF) 


This equation was developed in 1935 by C.F. Von Weizacker based on similarity of the nucleus with liquid drop 
as we mentioned earlier, and this equation is also known as Bethe-Weizacker equation[1]. It is based on 
consideration that total binding energy of the nucleus is a set of terms, each of term expresses one of 
characteristics of the nucleus, as follows[1] 


EB = YE, (4) 
i=l 


The terms of previous equation are given as following [4] 


- The first term is called volume energy term; it was named by this, because size of the nucleus is proportional 
to number of its nucleons. The interpretation of this term was based on existence of a specific average binding 
energy per nucleon for all nuclei, and it written as follows 


E, =, A (5) 
Where Ey is volume energy term and A, volumetric proportionality coefficient. 


- The second term is called surface energy term. This term was found based on phenomenon of surface tension 
in liquids. In a liquid drop, the molecules on surface of drop are weaker than the molecules inside it, and 
therefore their energy is greater. The same thing is achieved in the nucleus, where the nucleons located on 
surface is weaker than the nucleons located deep in the nucleus, and therefore it possesses additional surface 
energy, because inner nucleon exchanges influence with a greater number of nucleons than in state of outer 
nucleon located on surface of the nucleus, and this term is written as 


E, =-a, A” (6) 
Where E: is surface energy term and 4, surface proportionality coefficient. 


This term is added with a negative sign to binding energy equation because it decreases value of total binding 
energy. 
- The third term is Coulomb energy term, where existence of Coulomb force between protons leads to a 


decrease in value of binding energy. Coulomb potential is proportional to Z(Z — 1) , then expression Coulomb 
energy is written as follows 
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Ec =a. = (7) 
c C 
A 1/3 


where Eg is surface energy term and Aç Coulombic proportionality coefficient. 


This term is also added with a negative sign to binding energy equation, because repulsion energy decreases 
value of total binding energy. 

The previous three terms are considered classical terms. 

- The fourth term is called asymmetric energy term, and this term result from difference between number of 
protons and neutrons in the nucleus, as this difference leads to a decrease in stability of nucleus and thus a 
decrease in value of its total binding energy, so this term is also added with a negative sign to binding energy 


equation, and written as follows 
2 2 
N -Z A-2Z 
Eas =s I Fi ) =d;s í F ) (8) 


where Big is asymmetric energy term and @,, proportionality coefficient of the asymmetric energy. 


- The fifth term is called pairing energy term. This term was added because binding energy does not depend on 
number of nucleons only, but also on even or odd number of nucleons. It was found that the most stable 
nucleus in nature are even-even nuclei, and the least stable are odd-odd nuclei, the middle state is odd nuclei 
(even-odd or odd-even) and these nuclei are called odd nuclei or odd-A nuclei. The number of these nuclei is 
approximately one hundred nuclei found in nature. 

Accordingly, Fermi introduced a correction factor added to binding energy, which takes the following values 


- When A is even, N and Z even, the term is 0 = +a,A™", thus increases binding energy, because the 
nucleons are paired with each other. 

-When A iseven N and Z odd, the termis ô = —a, y. , thus binding energy is reduced. 

- When A is odd, this term does not take any value ô =O, and it is the middle state between two previous 


cases, as this term does not have any contribution to binding energy. This term can be written in abbreviated 
form as follows 


-3/4 
5 =(4,0)a, A (9) 
Where © is pairing energy term and a, coefficient of proportionality of pairing energy. An illustration of 


previous terms is shown in Figure (1). So, SEMF, takes the following form 


2 E 2 
e T (10) 


One of the most important methods used to calculate values of coefficients in SEMF is Least Square Fitting 


Ce a.o aO 
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Volume Surface Coulomb Asymmetry Pairing 


Source: https://commons.wikimedia.org/w/index php?title=File:Liquid_drop_model.svg&oldid=74367824. 


Figure (1): illustration for terms of SEMF in liquid drop model. 


5. LEAST SQUARES METHOD 


The least squares method is a standard method for finding solution to a set of equations that has more than one 
variable. This allows us to use this method to determine values of coefficients in theoretical formulas by fitting 
it with experimental values, where fitting is best possible when sum of squares of differences between 


experimental and theoretical values is smaller, as follows[2] 
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2 2 
=A TE) (11) 
Where f (x, )= p2 a, X A is a polynomial and 4, is coefficients of polynomial. 


Since function Y is dependent on several variables, it reaches a local maximum value at the point where all its 
partial derivatives of first order is zero, that is, when following equations are satisfied 


2 
Ox" g E 9, E 9, E0 (12) 


0a, Oa, Oa, ” a, 
The set of equations, whose number of variables is equal to number of its equations, is arranged in following 
matrix form 


(13) 


Bx Dap Exp = Ex {lal [Lev 
Uap Lap Dat. pxlt|fae]=|ZGFy) 


Elzy) equations, we get values of 
Ap sy Ay 5.,4,, and by 


n a xP á = Xx |. a 


By solving these T 2 g a . z 
pe pat g pa 


coefficients 


substituting them in polynomial f (x i ), we get required polynomial. 


6. ISOSPIN 


Heisenberg proposed this idea in 1932, when he considered protons and neutrons to be identical particles, he 
called nucleons, and assumed the very small difference between two masses resulted from Coulomb interaction 
between protons. Therefore, neutrons and protons represent two different states of one thing, which called 
nucleon, and this space is called isospin, and this space is accompanied by isospin operator, which we denote 


by T and takes following eigen value [5] 


T +D (14) 
The algebra of operators in isospin space is developed as in case of spin, isospin operator takes following 


form[5] 


1 
2° (23) 


Where T is Pauli isospin matrix, and accordingly, isospin space (or charge space) for a nucleon has two 


dimensions, like spin space. The third projection of isospin takes only two values for nucleons, which is [6] 
1 1 
T,=—1,=+-— (16) 
2 2 


The special value of nucleon according to Dirac's representation is IT Tas the special value of proton is 


l 1 1 
2? Oy 


—,—), while the special value of neutron is 

22 
The idea of isospin can be generalized for case of a nucleus that has Z protons and N neutrons, where it has a 
mass number of A=Z+N. In this case, the third component of isospin is added algebraically to each nucleon (for 


protons and neutrons) to get this compound for nucleus, as follows [6] 


T E An *(N Zj= (A 2Z ) (17) 


2 2 2 


As it is known, value of compounds cannot exceed the value of its vector, so value of isospin vector must be 
as[6] 


cS (18) 


N-Z 
2 
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The expectation value of two interacting nucleons is calculated from defining equation for total isospin which is 
equal to vector resultant isospin of each of two interacting nucleons as follows 


T =t, +t, (19) 
where f,,f,, is isospin of two interacting nucleons, by squaring both sides of previous equation, we get 


T? =t; +t,’ +2t,t, (20) 
1 
Given t, =t, = 5 . So, the expected value ff, is 


me De Ge) E (21) 


And since isospin is conserved for strong interactions (which corresponds to fact that the nuclear forces are 
symmetrical), it is therefore possible to adding a term related to isospin within total binding energy equation 


form[6], and thus term that can be added to SEMF given by equation (10) takes the following 
(c,t,)=2 T +1)-3 (22) 


7. SINGLE PARTICLE MODEL 


This model assumes that the nucleons in the nucleus of a specific type (protons or neutrons) are paired in a 
basic state with each other in form of pairs to give a total angular moment and a magnetic moment equal to 
zero, and these paired nucleons form a closed core. According to this, the total angular moment and isospin of 
even-even nuclei are equal to zero, and this is what we observe in nature, while the properties of odd-A nuclei 


returns to properties of unpaired nucleon[4]. 


8. CALCULATIONS AND RESULTS 


Calculations were made for 100 individual stable nuclei (odd-even and even-odd) located in different area 
(light-medium-heavy), and they are most of stable nuclei found in nature. 


8.1. values of coefficients in SEMF 
Depending on method of least squares fitting, where fitting is the best possible for binding energy equation 


given by equation (10) when sum of squares for differences between experimental and theoretical values is 
smallest, as follows 


0 


SL IDA (y, -EB (Z,,A,))=0 


Oa 
by Z.(Z, -1) (23) 
TA = 2D is —EB,,(Z,,A,))=0 


2 
oy Egy, -EB,,(Z,,A,))=0 


l 


Substituting equation (23) into (10), we get the following matrix equation 


ay cı | [bu Dy By bi 


as = c, | |B, ba Do, ba (24) 
Ac C3 || bz Dy, B33 ba 
ais C4 | [Par Ba Ba Pas 


Where values of elements [b ] are given as follows 
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bi =$ A; , bo =b,, =- A;", b =b; =-9 Z, (Z, =DA 
bq =b = -J 4, -2Z, : bn = DA Px =b, = Èz, (Z; =A 


(A, -2Z,)° Z;(Z,-1)° (25) 
by, =by = Dar „b3 = D a 


L 


Z,(A, -2Z,)° (A, -2Z,)* 
by =b; = m hu > a i 


i L i 


l 


The values of elements [c ] are given as follows 


Z (Z, -—1) (A. -2Z. 7 

— pae 2/3 — L L _ L L 

c= QA, C=) YA; i=- A pi A (26) 
By solving matrix equation (24), after substituting values of elements given in (25) and (26) with helping of a 
code designed for this purpose using famous program MATLAB, we get values of coefficients, which are given 


in units of MeV as follows 
a, = 15.8838, a, =18.1220, a, =0.7279, a,, =23.2012 


which is a new set of parameters for odd-A nuclei. 


8.2. values of coefficients in SEMF after adding isospin 
After adding isospin term given by (22), we add to set of equations (23) the following equation 


a (I )(y, -EB,,(Z,,A,))=0 (27) 
Thus, matrix equation (24) becomes as 
a | Le; | (by Bia ba Da Bic. t 
ds c, | | ba Ba by; ba bas 
ac |=] C3 |-| bz, By, B33 B54 Das (28) 
Aas C4 | [Pa Ba ba Bay Pas 
ap | Les | [Psi Ps: bss bsa bss | 


Where values of elements of matrix [b ] are given as follows 
2 5/3 2/3 
b =) A; db, =b, =—-9 A; » b =b; == ZZ, -DA; 


b =b; =); -2Z,} ‚b; =b, = 2A; [ar T +1)-3] 


2 
by, =} 47", b,, =b; =).2,Z, -DA b, =b, apan a 


- AY 
Z?(Z, -1° 
b»; =b; =- A" [rw +03]; ba = 2o A 
A EEA Z (Z, -0 
Ds, = bay = a bs =b; = 2 A [2T T +1)-3] 


Dy, =y = „bas =bs4 =-)) eee 2A [2T T +1)-3] 
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The values of elements of matrix [c] are given as follows 


Z, (Z, -1 A, -2Z, 
c=) yA, Cy =-)_ AP” as +) < ue Yy, í i 2 
i i i i i i (30) 
c5=-) [We +1)-3]y, 


By solving matrix equation (28), after substituting values of elements given in equations (29) and (30), we get 
values of coefficients a, ‚aș ‚ac ,d,5 ,4, , which are given MeV as follows 


a, =16.8290, a, =22.0231, a, =0.7770, ay; =25.0149, a, =8.3892 


We calculated standard deviation between our calculated binding energy and experimental values for on 
hundred stable odd-A nuclei, as we use the following standard deviation relation[7] 


1 Š i i] 
o= aks -E (31) 


B (theo .) B (exp.) | 


We found value of standard deviation for these nuclei without adding isospin term is © = 1.831(Mev ), and 


we found standard deviation value after adding isospin term is O = 1.662(Mev ) , that is, adding isospin term 


for SEMF reduces value of standard deviation. And it appears in Figure (2) distribution of standard deviation 
values of nuclei calculated by us before and after adding isospin in terms of mass number. 
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Figure 2: distribution of standard deviation values of our calculated nuclei before adding isospin 
(dotted line) and after adding it (solid line) in terms of mass number. 


When comparing standard deviation value for calculated nuclei with previous studies, we found value, which 
we calculated after adding isospin term, is better than values calculated in previous studies mentioned in 
references[8-10], which depending on different methods and different approximations, and it is better than 


values calculated in reference[11], which was performed by using coefficients deduced by Evans and Arya, 
using UAEN method, and these comparisons are shown in table (1). 
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Table (1). Comparison between standard deviation value calculated by us with previous studies 


References [10] [11] [11] [11] | Present Work | Present Work 
UAEN | Arya | Evans Without iso. With iso. 
o(MeV ) |4711 | 3.888 | 3.795 | 1.9307 | 3.070 | 2.729 | 2.412 1.831 į 1.662 


We also found that our standard deviation value is better than study mentioned in reference[12], Which 
depending on adding different terms to SEMF, which are, exchange Coulomb term, Wigener term, surface 
symmetry term, pairing term, curvature term and shell effect term, and also calculated with effect of adding all 
the previous terms. So, we notice adding isospin effect term is great importance compared to adding other 
terms. These comparisons are shown in table 2. 


Table 2. Comparison of our calculated standard deviation value with values of meee 


Exchange Surface | Pairing | Curvature Shell Present į Present 
Coulomb symmetry effects | terms Work Work 
Without 

iso. 


9. CONCLUSION 


Using a code designed using MATLAB program, we calculated values of four coefficients of SEMF, which were 
as follows 


a, =15.8838, a, =18.1220, a, =0.7279, ay, =23.2012 


for one hundred stable nuclei that have an odd mass number (even - odd and odd- even) nuclei, using least 
squares fitting method. We also added another term to SEMF related to isospin using single-particle model. And 
we found values of coefficients of SEMF after adding this new term, which were as follows 


a, = 16.8290, a, = 22.0231, a, =0.7770, a,, = 25.0149, a, =8.3892 


We found adding isospin term to SEMF reduces value of standard deviation between binding energies of our 
calculated nuclei and experimental values. 

When comparing values of binding energies for studied nuclei with values of previous studies, we found that 
our calculated values are better than values in these studies. 

The addition of isospin to SEMF can be tested for nuclei having an even- even nuclei, at which point calculations 
must be performed by adding pairing term, and a way to find out isospin of interacting nucleons. 
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